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Monopole field configurations have been extensively studied in both Abelian and non-Abelian 
gauge theories. The question of the quantum corrections to these systems is a difiicult one, since 
the classical monopoles have non-perturbatively large couplings, which makes the standard, pertur- 
bative methods for calculating quantum corrections suspect. Here we apply a modified version of 
Heisenberg's quantization technique for strongly interacting, nonlinear fields to a classical solution 
of the SU(3) Yang-Mills field equations. This classical solution is not monopole-like and has an 
energy density which diverges as r — > co. However, the quantized version of this solution has a 
, monopole-like far field, and a non-divergent energy density as r ^ oo. This may point to the con- 

' elusion that monopoles may arise not from quantizing classical monopole configurations, but from 

^S) ' quantizing field configurations which at the classical level do not appear monopole-like. 

o 

Q 

m 
(N 



> 
oo 



X 



I. INTRODUCTION 



Monopoles have been studied within the context of both Abelian ||] and non-Abelian Q gauge theories. Usually, the 
approach is to start with a classical monopole configuration {i.e. having a magnetic field which becomes Coulombic 
as r — » oo) and then consider the quantum corrections to the system. One difficulty with this approach is that the 
coupling strength of monopole configurations is large due to the Dirac quantization condition Q which requires that 
' there be an inverse relationship between electric and magnetic couplings. A perturbatively small electric coupling 
requires a non-perturbatively large magnetic coupling. The perturbative quantization techniques do not work well 
with monopoles for the same reason that they have trouble with QCD in the low energy regime : the couplings are 
I non-perturbatively large. 

^ ' Here we will use a modification of a non-perturbative quantization method originally used by Heisenberg to quantize 
^ , nonlinear, strongly coupled spinor fields ||] Q. However, rather than applying this quantization to an SU(3) field 
configuration which is monopole-like already at the classical level, we apply it to a non-monopole configuration, 
Oh which has fields and an energy density which diverge at spacial infinity (r — > oo). After applying Heisenberg's 
^ ' quantization technique to this solution we find that the fields and energy density become physically well behaved, 
and the asymptotic magnetic field becomes monopole-like. This result may indicate that monopoles are inherently 
quantum objects. Rather than quantizing a classical field theory configuration which appears monopole-like, real 
monopoles might arise as a consequence of quantizing non-monopole classical solutions. 



II. CLASSICAL SU(3) SOLUTION 

In this section we will briefiy review the classical, spherically symmetric, SU(3) Yang-Mills theory solution to 
which we will apply the modified Heisenberg quantization method. We begin with the following ansatz for pure SU(3) 
Yang-Mills theory 

A, ^ iA'^(A^,.+A,",)^u;(r), (1) 
Af = (A^,- - A,",) + A,", {eu,x' + ea,x^) -3 .(r), (2) 



here A" are the Gell - Mann matrices; a — 1,2, .. . , 8 is a color index; the Latin indices i,j,k,l = 1,2,3 are space 
indices; ~ —1: r, 9, (p are the usual spherically coordinates. This is a simplified version of an ansatz considered by 
several groups 
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Under this ansatz the Yang - Mills equations {D^F°- 1^^^ = 0) become 



r^v" 



r^w" 



' V — vw , 



6wv^ . 



(3) 
(4) 



where the primes indicate differentiation with respect to r. In the asymptotic limit r 
approach the form 

V w Asm{x°' + 0o) I 

„ a - 1 cos {2x"' + 2(j)o) 



oo the solutions to Eqs. 



w « ± 

i2 



= a(a-l) 



(5) 
(6) 
(7) 

where x = r / ro is a dimensionless radius and rp, 0O: and A are constants. When Eqs. were studied numerically 

it was found that the approximate form of the solution in Eqs. (||) - (|7|) was a good representation even for x fairly 
close to the origin. It was also found that for a broad range of boundary conditions for the functions v, w that a tend 
to fall in the range 1.2 < a < 1.8. Thus for large x the ansatz function w is a smoothly, increasing function, which 
yields an increasing time-like gauge potential from Eq. (|l|). This can be compared with the confining potentials used 
in some phenomenological studies of quarkonia bound states It is also possible to find analytic solutions with 
such increasing gauge potentials if one couples the Yang-Mills fields to a scalar field jl^. However, this increasing of 
w also leads to the field energy and action of this solution diverging as a; — > oo. 

The "magnetic" (-fff = Fjf^) and "electric" {Ef = Fq^) fields associated with this solution can be found from the 
non-Abelian gauge potentials, A^, and have the following behavior 



oc 



rw' — w 



oc w', 
E^ oc — , 



^ r r 

here for E^, H^, and the color index a = 1, 3, 4, 6, 8 and for H^, E^ and E^ a = 2, 5, 7. 



(8) 
(9) 



III. HEISENBERG'S NON-PERTURBATIVE QUANTIZATION 

Although the confining behavior of the above classical solution is of interest due to its similarity with certain 
phenomenological potentials, its infinite field energy makes its physical importance/meaning uncertain. One possible 
resolution to this classical solution's divergent field energy is if quantum effects removed its bad long distance behavior. 
The difficulty is that strongly interacting, nonlinear theories are notoriously hard to quantize. In order to take into 
account the quantum effects on this solution we will employ a variation of the method used by Heisenberg |^ in 
attempts to quantize the nonlinear Dirac equation. We will outline the key points of the method by using the 
nonlinear Dirac equation as an illustrative example. The nonlinear spinor field equation considered by Heisenberg 
had the following form : 

j'^d^.i^ix) - l^^iPiU)] = (10) 

where 7'^ are Dirac matrices; 'ijj{x), ip are the spinor field and its adjoint respectively; ^[iplipip)] is the general nonlinear 
spinor self interaction term which involved three spinor fields and various combinations of 7^'s and/or 7^'s. The 
constant / has units of length, and sets the scale for the strength of the interaction. Next one defines r functions as 

T(xiX2...|2/i2/2...) = {0\T[^ixi)^ix2)...riyiW{y2)-m (11) 

where T is the time ordering operator; |$) is a state for the system described by Eq. (p^. Applying Eq. ( pX| ) to ( pl| ) 
we obtain the following infinite system of equations for various r's 

d 

^~'^^ir)jn^ T(a;i---a;„|yi...y„) = ^[T{xi...XnXr\yi..-ynyr)] + 

5{Xr - yi)T{xi...Xr-lXr+l--Xn\y2--yr-iyr+l--yn) + 

S{Xr - y2)T{xi...Xr^lXr+l...Xn\yiy2--yr-iyr+l--yn) + •■• (12) 
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Eq. ( p^ represents one of an infinite set of coupled equations which relate various order (given by the index n) of the 
T functions to one another. To make some head way toward solving this infinite set of equations Heisenberg employed 
the Tamm-Dankofl[ method whereby he only considered r functions up to a certain order. This effectively turned the 
infinite set of coupled equations into a finite set of coupled equations. 

For the SU(3) Yang-Mills theory this idea leads to the following Yang-Mills equations for the quantized SU(3) gauge 
field 







(13) 



here F"'^^'^ is the field operator of the SU(3) gauge field. 

One can show that Heisenberg's method is equivalent to the Dyson-Schwinger system of equations for small coupling 
constants. One can also make a comparison between the Heisenberg method and the standard Feynman diagram 
technique. With the Feynman diagram method quantum corrections to physical quantities are given in terms of 
an infinite number of higher order, loop diagrams. In practice one takes only a finite number of diagrams into 
account when calculating the quantum correction to some physical quantity. This standard diagrammatic method 
requires a small expansion parameter (the coupling constant), and thus does not work for strongly coupled theories. 
The Heisenberg method was intended for strongly coupled, nonlinear theories, and we will apply a variation of this 
method to the classical solution discussed in the last section. 

We will consider a variation of Heisenberg's quantization method for the present non-Abelian equations by making 
the following assumptions [ pi] : 

1. The physical degrees of freedom relevant for studying the above classical solution are given entirely by the two 
ansatz functions v, w appearing in Eqs. (^), (^. No other degrees of freedom will arise through the quantization 
process. 

2. From Eqs. (^, (^ we see that one function w{r) is a smoothly varying function for large r, while another 
function, f (r), is strongly oscillating. Thus we take w{r) to be an almost classical degree of freedom while v{r) 
is treated as a fully quantum mechanical degree of freedom. Naively one might expect that only the behavior of 
second function would change while first function stayed the same. However since both functions are interrelated 
through the nonlinear nature of the field equations we find that both functions are modified. 

To begin we replace the ansatz functions by operators w{x)^ v{x). 



r'v" 



V — vw 



These equations can be seen as an approximation of the quantized SU(3) Yang-Mills field equations (n3| 
into account assumption (2) we let w ^ w become just a classical function again, and replace t)^ in Eq, 
expectation value to arrive at 



(14) 
(15) 

Taking 
(|15|) by its 



2 -// 

r V 



vw 



where the expectation value 
Eq. (|l|) to get 



6w{v ) 

is taken with respect to some quantum state \q): {v'^) — {q\v'^\q). 



(16) 
(17) 

We can average 

(18) 
(19) 

Eqs. ( p^ ( |l9| ) are almost a closed system for determining (v) except for the {v'^) and (v^) terms. One can obtain 
differential equations for these expectation values by applying r'^d/dr to v"^ or ■0^ and using Eqs. (16) - (p^). However 



= - iv) {v)i 



the differential equations for (w^) or {v'^) would involve yet higher powers of i) thus generating an infinite number of 
coupled differential equations for the various (f"). In the next section we will use a path integral inspired method 
p2[ to cut this progression off at some finite number of differential equations. 



IV. PATH INTEGRATION OVER CLASSICAL SOLUTIONS 



Within the path integral method the expectation value of some field $ is given by 



(20) 
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The classical solutions, <I>c/, give the dominate contribution to the path integral. For a single classical solution one 
can approximate the path integral as 

ye^S[<[.]^^^^giS[4.<=,] (21) 

where A is a normalization constant. Consequently the expectation of the field can be approximated by 

$e^s[<E.]^^ « ^^i. (22) 

We are interested in the case where $ is the gauge potential j which case our approximation becomes 



M^>« / Ul)^ e ^'-olofA;) (23) 

(v4°)^Q are the classical solutions of the Yang - Mills equations labeled by a parameter ^o- In the present case the 
classical solution with the asymptotic form (||) ^ has an infinite energy and action. When one considers the Euclidean 
version of the path integral above the exponential factor in ( |2^ ) becomes exp[— S'[(^°)0(,] which for an infinite action 
would naively imply that this classical configuration would not contribute to the path integral at all. However, 
there examples where infinite action classical solutions have been hypothesized to play a significant role in the path 
integral. The most well known example of this is the meron solution [T^ , which has an infinite action. Analytically 
the singularities of the meron solutions can be dealt with by replacing the regions that contain the singularities by 
instanton solutions. Since instanton solutions have finite action this patched together solution of meron plus instanton 
has finite action. However, the Yang-Mills field equations are not satisfied at the boundary where the meron and 
instanton solutions are sewn together. In addition recent lattice studies |Q have indicated that merons (or the 
patched meron/instanton) do play a role in the path integral. Here we will treat this divergence in the action in an 
approximate way through a redefinition of the path integral integration measure. Since the divergence in the action 
comes from first term of the w{r) ansatz function of (|^), which does not contain 4>q we will take the action for different 
's to be approximately the same : S[{A'^^)^g\ Sq oo. Then changing the functional integration measure in the 
following way ~ D{A'^)^g e~''^°I?(A° )0q - allows us to approximate the expectation value of as 

E {^"Xp^o (24) 

\ ■'00 



over all 
classical solutions 



where is the probability for a given classical solution. We will consider the classical solutions whose asymptotic 
behavior is given by (^)-(|^), and we will take the different solutions (as distinguished by different 0o's) to have equal 
probability p^g w const. Therefore 

(y) ~ L f y^^d(bo I sin(x" + 0o)#o = 0, (25) 



27r J " 27r 
b 

l^y^) ^ Lj yld<jy^ ^^J sin2(a;" + 0o)rf0o = ^, (26) 


2tt 2-n 

{v') ~ ^cid(bo = ^ y sin3(a;" + </>o)d0o = 0, (27) 


2tt 2ir 



(„4^ « i_ 1 „4 ^ ^ 1 gi^4(^a ^ ^ 3^4 (28) 



here Vd is the function from the Eq.(||). The path integral inspired Eqns. (25) - (|28| ) are the heart of the cutoff 
procedure that we wish to apply to Eqns. (p^ , (|T9|). On substituting Eqs. (25), (|26|) into Eqs. (|T^), (|9|) we find 
that Eq.(|l^) is satisfied identically and Eq.(|19|) takes the form 

r^w" = 3A^w = a{a - l)w, a > 1 (29) 
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which has the solutions 

w = War", (30) 

^ = ^ (31) 

where wq is some constant. The first solution is simply the classically averaged singular solution (^ which still has 
the bad asymptotic divergence of the fields and energy density. The second solution, (31), is more physically relevant 
since it leads to asymptotic fields which are well behaved. 

The solution of Eq. ( plj ) implies the following important result: the quantum fluctuations of the strongly 
oscillating, nonlinear fields leads to an improvement of the bad asymptotic behavior of these nonlinear 
fields. This means that after quantization the monotonically growing and strongly oscillating components of the 
gauge potential become functions with good asymptotic behavior. 

As r — *■ oo wc find the following SU(3) color fields 



— « — with <3 = g' 

(i?^.,) « {v')^0, (33) 



(H^) oc «^ with Q = -a{a-l)-l, (32) 



(E;) „ . (34, 

« ^ = »^ (35) 

We can see that as r ^ oo K-E^)! ^ particular at infinity we find only a monopole "magnetic" field 

Hr ~ Q/f'^ with a "magnetic" charge Q. This result can be summarized as: the approximate quantization 
of the SU(3) gauge field (by averaging over the classical singular solutions) gives a monopole-like 
configuration from an initial classical configuration which was not monopole-like. We will call this a 
"quantum monopole" to distinguish it from field configurations which are monopole-like already in the classical 
theory. 



ENERGY DENSITY 



The divergence of the fields of the classical solution given by Eqs. (||) - (0) leads to a diverging energy density for 
the solution, and thus an infinite total field energy. The energy density e of the quantized solution is 



The first two terms on the right hand side of Eq. (|36|), which involve the "classical" ansatz function g{r), go to zero 
faster than l/r'* as r ^ oo due to the form of g{r) in Eq. (^Tj). Thus the leading behavior of e is given by the last 
two terms in Eq. (^) which have only the "quantum" ansatz function f{r). To calculate (/'^) let us consider 

r'^ {f'{r')f'{r)} = {f'{r')fir)) (/'(r')/(r)) (l + g'ir)) (37) 

In the limit r' ^ r we have 

r' {f'{r)f"{r)) = ^ {f'\r))' = \ {f\r))' \ (f^r))' (l + g^r)) . (38) 



From Eqs. ( p6D and ( p8| ) we see that (/' )' = which implies (/' } = const. Thus the third term in ( p6D gives the 
leading asymptotic behavior as r — > oo to be 

e « ^ (39) 

and the total energy of this "quantum monopole" (excitation) is infinite. This fact indicates that our approximation 
( pO| ) is good only for (/"(r)) calculations but not for the derivative (/'^(t")). . 
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VI. CONCLUSIONS 



Starting from an infinite energy, classical solution to the SU(3) Yang-Mills field equations we found that the bad 
asymptotic behavior of this solution was favorably modified by a variation of the quantization method proposed by 
Heisenberg to deal with strongly coupled, nonlinear field theories. In addition, although the original classical solution 
was not monopole-like, it was found that the quantized solution was monopole-like. This may imply that if real 
monopoles exist they may be inherently quantum mechanical objects i.e. that monopoles arise from the quantization 
of non-monopole classical solutions. This is to be contrasted with the standard idea that monopoles result from 
quantizing solutions which are already monopole-like at the classical level. One possible application of this is to the 
dual-superconductor picture of the QCD vacuum. In this picture one models the QCD vacuum as a stochastic gas of 
appearing/disappearing monopoles and antimonopoles as in Figj^. These monopole/antimonopole fiuctuations can 
form pairs (analogous to Cooper pairs in real superconductors) which can Bose condense leading to a dual Meissner 
effect []l6| jl7| expelling color electric flux from the QCD vacuum, expect in narrow flux tubes which connect and 
confine the quarks. Lattice calculations confirm such a model [ p^ : monopoles appear to play a major role in the 
QCD lattice gauge path integral. Based on the results of the present paper it may be that the monopoles which 
are considered in the dual superconductor QCD vacuum picture should be the "quantum" monopoles discussed here 
rather than "classical" monopoles {i.e. monopoles which are already monopole-like at the classical level). 



quantum monopoles/antimonopoles 




FIG. 1: QCD vacuum ~ stochastic gas of quantum monopoles/antimonopoles. 
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